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Abstract 

For Bloch electrons in a magnetic field, explicit solutions are obtained at the 

center of the spectrum for the Bethe ansatz equations recently proposed by 

Wiegmann and Zabrodin. When the magnetic flux per plaquette is 1/Q where 

Q is an odd integer, distribution of the roots is uniform on the unit circle in 

the complex plane. For the semi-classical limit, Q — > oo, the wavefunction 

obeys the power low and is given by |?/>(£)| 2 = (2/ sin-7rx) which is critical and 

unnormalizable. For the golden mean flux, the distribution of roots has the 

exact self-similarity and the distribution function is nowhere differentiable. 

The corresponding wavefunction also shows a clear self-similar structure. 
PACS numbers: 72.15.Gd, 03.65.Fd 
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The Azbel-Hafstadter-Wannier problem of two dimensional Bloch electron in a magnetic 
field is an old problem that often brings new excitements. The equations of motion can be 
reduced to one-dimensional ones, which appear in many different physical contexts, ranging 
from the quantum Hall effect |L| to quasi-periodic systems (for a review, see for example 
ref. 0.) Its topological character has been revealed both for periodic boundary conditions 
and for systems with edges ||. The interplay of two intrinsic periods, the period of the 
lattice and that of the magnetic flux is essential in this problem. Previously many studies 
were done for the cases when the two periods are commensurable, i.e. the magnetic flux 
per plaquette, 0, is rational. But the incommensurate cases (with an irrational 0) are more 
interesting in that the spectrum is known to have an extremely rich structure like the Cantor 
set and to exhibit a multifractal behaviors. |4],|2| Another interesting case is the weak field 
limit: When the flux is small, the semi-classical treatment of the WKB type is justified. 
Although some properties of the wavefunctions are known [||||, much of the knowledge we 
have came from numerical studies for rational fluxes. [[| In this letter, we derive several 
analytical results for the weak-field limit and for the incommensurate golden-mean flux, by 
following the highly innovative Bethe ansatz approach recently proposed by Wiegmann and 
Zabrodin [|7| using quantum group techniques. 

The Hamiltonian is given by H = J2 m ,n( c ln+i,n e ' 9 " 1, ° c m,n + c ln,n+i eld " l ' nc m,n + H.c.) where 
c mtn is the annihilation operator for an electron at site (m,n). We first assume the flux 
per plaquette is rational: = P/Q with coprime odd integers P and Q, and later take 
appropriate limits of 0. Let us take the diagonal gauge, 0^ nn = +7r0(n+m), Q v m n = — 7r0(n + 
m). || Then the Schrodinger equation for the one particle state \^f) = J2 m ,n Vv^m n|0) 
is written as i(q l+1 + g^^+i — i^q 1 ^ 1 + q~ l )tyi_i = E^i at the mid-band point where 
"0m,n = i m+n ^> m+n and q = e tn( ^ = e l7rP/,< 3. Wiegmann and Zabrodin observed J7[ that if we 
write = ^(q 1 ) with ^(z) a polynomial of degree Q — 1, this difference equation coincides 
with the eigen-equation of a combination of generators of the quantum group U q (sl2) in the 
representation provided by ^{z). Thus the zeros, {zi}, of the polynomial ^{z) must satisfy 
the Bethe ansatz equations (BAE) |7j 
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qzf + 1 ^ zi - qz m 

For the zero energy E — 0, they showed that ^(z) is given by the so-called contin- 
uous q-ultraspherical polynomial |Ll]] as ^(z) = j^y-{—iz) n P n (—iz) , where P n (-2) = 
Efe=o ( g;g:j;$);! fc ^- 2fc , n = (Q - l)/2 and (a; g) fc = Efci(l - <). 

In order to understand the properties of the wavefunctions at E = 0, the center of the 
spectrum, we try to solve the BAE ([I]) explicitly The polynomial P n {z) satisfies a difference 
equation ]12| 



(1 - qz 2 )P n (qz) + (q- z 2 )P n {q- x z) = 0. (2) 

Put z = q ±x l 2 1 then we get P n (q ±3 ^ 2 ) = 0. So g ±3 / 2 are roots. By iteration, the complete set 
of roots of the BAE are given by 

z m = iq 2m ~ Xl \ iq- 2m+l ' 2 , m = 1, • • • , (Q - l)/2. (3) 

They are all on the unit circle. Let us write the roots as z m = e l9m and consider the 
distribution of 9 m . 

Consider first the special case P = 1. The roots {z m } distribute uniformly on the unit 
circle except near z = ±i. The roots for Q = 21 are shown in Fig. [TJ. In the semi-classical 
limit Q — > oo, that is, q — > 1, the distribution function p(8) = limg^oc QA6 is smooth 
(constant) where A8 = 9 m+ i — 9 m . A continuous behavior of p(8) is usually obtained in the 
exactly solvable models (the Heisenberg chain, the Hubbard chain, etc.) in which p(8) is 
determined by an integral equation. 

When the flux is irrational, the situation is quite different. Take the flux = 1/r = 
(V5 - l)/2, where r is the golden mean. To reach this flux, we consider a sequence of 
rational fluxes <pk — Pk/Qk-, where Qk = F 3 ^ + i, P^ = F^ and is a Fibonacci number 
defined by Fk+i = Fk + F^-i, F\ — 1, F — 1. In this case, the two types of roots in 
Eq. are nested. To gain an insight into the distribution of roots, it is helpful to consider 
the pseudo roots which are defined also by Eq. (^) but with the range of m modified to 



m = —{Q — l)/2, • • • , 0. In Figs. |2] we show the distributions of the roots (black points) 
and pseudo roots (gray points) for several ( k = 1,2,3,4,5,6). Here the radius of the 
unit circle has been scaled so as to show all the cases at once. These figures clearly show 
that there is a branching rule for the true roots (denoted by A) and pseudo roots (by B) as 
follows: 

A 3 -> A 3 B 2 A 3 B 2 A 3 

A 2 -> A 3 B 2 A 3 

B 3 -> B 3 A 2 B 3 A 2 B 3 

B 2 -> B 3 A 2 B 3 . (4) 

The initial condition is B 3 A 2 B 3 A 2 (cyclic). In the fc-th stage of the sequence, the number of 
clusters of the true roots A 3 and A 2 are Qh-x — 1 and Ph-x + 1 respectively. This branching 
rule gives rise to a self-similar structure for the distribution p(6) in the limit k — > oo. 
To characterize the distribution, let us define the generation of a root. According to the 
branching rule (f|), each true (pseudo) root branches into a cluster of 3 new true (pseudo) 
roots, each of which in a sense has a parent. At the same time, between these clusters of new 
pseudo (true) roots, there is a pair of new-born true (pseudo) roots which have no parent. 
We assign the generation number to a root so that it is 1 when the root does not have a 
parent, otherwise it is one plus that of its parent. Let us denote the number of true (pseudo) 
roots in the k-th stage with generation g by nA{g, k) (riB{g,k)) (g — 1, • • • , k). (Then in 
the special case P = 1, n^g^k) = ns{g,k) = 8k g -) In the present case, we get a recursion 
formula by the branching rule as riA{g, k) = 3riA{g — 1, k — 1), g = 2, • • • , k, n^(l, k) = 
2(P fc _x + 1). Thus n A (g, k) = 2 ■ ^{Pk-g + 1) and n B (g, k) = 2 • &- x {P k _ g - 1). The above 
considerations exemplify the difference of the distributions between the semi-classical limit 
and the incommensurate case. The distribution of the roots has a self-similar structure and 
the function p{9) is nowhere differentiable in the incommensurate limit, while that for the 
semi-classical limit is smooth. We believe this is characteristic to the incommensurate case. 
Another way to characterize the distribution is to map to the dual (reciprocal) space. 
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This can be done for arbitrary P and Q. We lift the 9 m to the real axis periodically. On 
the real axis , the true and pseudo roots occupy a lattice of points {j/2Q | j: integer} with 
spacing 1/2Q. Thus we perform the Fourier transform by Sq(p) = J^JL-oo eip SqU) where 
Sq(J) is the so-called defining function: Sq(J) = 1 if there is a true root at j/2Q, otherwise 
S Q (j) = 0. Then 

S Q (p) = ^j: s Q 5 (P-Pr), (0<p<2vr), (5) 

where s° Q = (Q - 1), s r Q = (-) r+1 [cos(iPp r )]- 1 (r = 1,---Q - 1), and p r = 2nr/Q. 
In the semi-classical limit, \Sq(p)\ 2 is well defined and behaves smoothly. On the other 
hand, \Sq(p)\ 2 is not even differentiate in the incommensurate limit since P — > oo. Also 
it can be shown that the original defining function is given by Sq(j) = (^)Xm=i[l — 
(-) n cos(27rjVQ)/cos(7rPn/Q)]. 

Next let us consider the wavefunction. Using the explicit roots we found above, the wave 
function at site j can be written in a compact factorized form 

% = ( <D >r«! q 2 ) iQ - m ^- 3/2 ; q-\ Q -m- (6) 

It is convenient to shift the site, j = j + jo, by an amount jo where jo is determined by 
Pj = {Q-p) /2 (mod 2Q ). Then % = at ] = 2m, -2m + 1 (m = 1, • • • , (Q - l)/2), and 
the wavefunction is nonzero only at j — 1, 3, - • • , Q and Q + 1, Q + 3, • • • 2Q. Using Eq. (||), 
the amplitude of the wavefunction is calculated as 

(Q; rV /2 u . ir(2m-])P . ix(2m + j - 1)P |2 . . 

First let us present the exact results explicitly for the special case P = 1. A direct and 
careful calculation leads to a very simple result for the amplitude 

l^)| 2 = — (0<x<l), (8) 
sm(7rx) 

up to a constant factor, where if)(x) = ^j, x — (2j — 1)/2Q, and the semi-classical limit 
Q — > oo is taken. Thus the squared amplitude of the wavefunction is given by the inverse 



chord distance in the semi-classical limit. The recursion relation \^]\ 2 = |\l/j_2| 2 sin 2 (7r(j — 
2)/Q)/ sin 2 (7r(j — l)/Q) obtained from Eq. (|7|), has played a key role. [13| For a finite Q, a 
correction factor appears only near the edges (x = and 1). For example, near the edge x ~ 
0, the finite size correction is given by \ip(x2i+i)\ 2 = C(Z)(2/ sin7r£2Z+i), with C(l) = tt(1 + 
V^nLiU - l/ 2 *0 2 > wh ere C(0) = tt/4 = 0.78539..., C(l) = 5tt/16 = 0.98174..., C(2) = 
8l7r/256 = 0.99402..., •• •. So the finite-size correction factor C(Z) converges to unity very 
rapidly and the Eq. (H) is quite accurate even at small Q. The norm of the wavefunction 
is logQ + const, and unnormalizable which is characteristic to a critical wavefunction. In 
Fig. |3], the amplitudes of the analytic wavef unctions, normalized by the peak height, are 
shown for several values of Q. 

Next let us discuss the case with golden-mean flux. We plot the analytic results Eq. ([?]) 
in Fig. f| for a sequence of rational fluxes converging to 1/r. One can easily recognize 
the self-similar behavior of the wave function. Each peak branches into three peaks in the 
next stage. Presumably these are the reflection of the self-similar distribution of the roots, 
Eq. (^). The multifractal analysis [Q] done numerically gives a smooth f(a) shown in Fig. 
This clearly shows that this wavefunction is multifractal and critical. We note the striking 
resemblance of these wavefunctions to that of the Id quasicrystal Fibonacci lattice at the 
center of the spectrum. fl4| An analytical derivation of f(a) is under progress. The latter 



was obtained exactly by a different technique and f{a) is obtained analytically. 

In conclusion, we have found analytical solutions to the BAE ([!]) that describes the Bloch 
state in a magnetic field with zero energy. The flux per plaquette is = Pj Q with coprime 
odd integers P and Q. All roots are on the unit circle, and the defining function for the 
roots is explicitly derived. 

When P = 1, the roots distribute uniformly for any odd Q. In the semi-classical limit 
Q —>■ oo, the density function of roots is smooth. When the flux (p are ratios of the successive 
Fibonacci numbers which converge to the golden mean, we found the branching rule for the 
roots, which makes the density function p of the roots exactly self-similar and nowhere 
differentiable. 



We are able to get explicit wavefunctions too. In the semi-classical compact 
expression is derived for the critical wavefunction, which turns out not normalizable. The 
squared amplitude, ^(x)! 2 , is given by the inverse chord distance, 2/ sin(7ra;). For the golden 
mean flux, the wavefunction also has a clear self-similar branching structure. 

Finally we remark that the solutions (^) to the BAE has a relatively simple structure, 
and positions of the roots are not affected by a finite Q correction. These facts together 
with the simple form (H) of the amplitude of the wavefunctions suggest that the symmetry 
of the problem, at least at the center of the spectrum, is higher than originally expected. 
It might just be the quantum group Uq^slz), though a clear-cut proof is still lacking. The 
situation would be similar to the Haldane-Shastry model where the symmetry is known to 
be higher than what is naively expected. |l6| . [T7f 
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FIGURES 



FIG. 1. Roots of the Bethe ansatz equation for the special case P = 1 and Q = 21. 

FIG. 2. Roots of the Bethe ansatz equation for the ratio of the Fibonacci numbers. 
<f>k = Pk/Qk = 3/5,13/21,55/89,233/377,987/1597,4181/6765 (a) in the whole complex plane and 
(b) enlarged figure. In each case, the roots are always on the unit circle. We scaled the radii to 
show the branching rule clearly. 

FIG. 3. Squared amplitudes for the wavefunctions: the special cases P = 1 and 
Q = 5, 21, 89, 377. The wavefunctions are normalized by the peak heights. 

FIG. 4. Amplitude of the wavefunctions: for the ratios of successive Fibonacci numbers: 
(a) whole region ((f) = 3/5, 13/21, 55/89 ) and (b) enlarged plot in the region near 0.89. 
((f) = P k /Q k = 3/5, 13/21, 55/89, 233/377, 987/1597, 4181/6765) The wavefunction with larger 
value of Q is shaded darker. The wavefunctions are normalized by the peak height. 

FIG. 5. The plot of f(a) for the wavefunction in the golden mean case. 



10 



0.5 



-0.5 



• 
• 


\ 

• 
• 

• 
• 
• 


: 
• 
• 
• 
• 
• 


• 
• 
• 
• 
• 
• 



-1 -0.5 

Fig.1 



0.5 1 



Fig.3 




JL 




oait 



Hi 



J 



(a) 



0.5 



-0.5 



-0.5 



Fig.2(a) (b) 

0.74 



0.735 
0.73 

0.725 
0.72 

0.715 
0.71 



0.675 

Fig.2(b) 



1 

0.8 
0.6 
0.4 
0.2 



0.85 0.86 0.87 0.88 0.89 0.9 0.91 0.92 

Fig.4(b) 

Y. Hatsugai, M. Kohmoto, and Y. S.Wu 



0.5 
Fig.5 



